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Abstract 

We study a boundary value problem associated with a system of two second order differential equations 
with cubic nonlinearity which model a film of superconductor material subjected to a tangential magnetic 
field. We show that for an appropriate range of parameters there are asymmetric solutions, and only 
trivial symmetric solutions. We then correct an error of the authors in [9] and show that the associated 
energy function is negative for the asymmetric solutions, and zero for the trivial symmetric solution. 
It follows that a global minimizer of the energy is asymmetric. This property resolves a conjecture of 
Marcus [13]. 
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Introduction 



In this paper we continue our recent studies ([9], [10]) of the one dimensional Ginzburg-Landau model 
[8] for superconductors. Our main objective is to investigate the model for the existence of asymmetric 
minimizers of the appropriate energy integral. It is expected that the physically interesting solutions will be 
energy minimizers. These solutions satisfy a symmetric boundary value problem which was known to have a 
set of symmetric solutions. We will show, however, that in one space dimension, for some parameter values 
the energy minimizer is an asymmetric solution to this symmetric boundary value problem. 

The problem may be compared to studies in two dimensions, where the symmetry of the energy minimizer 
is unresolved. Recent papers on this include [12] and [14]. Work in two dimensions has largely dealt with a 
reduced problem in which variations in the magnetic field are ignored. In contrast, it is possible to include 
the magnetic field as a variable in the analysis of the one dimensional model. 

In order to properly describe our results we first need to give a brief development of the problem as well 
as a summary of our previous investigations. 
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In 1950 Ginzburg and Landau [8] proposed a model for the electromagnetic properties of a film of 
superconducting material of width d subjected to a tangential external magnetic field. Under the assumption 
that all quantities are functions only of the transverse coordinate, they proposed that the electromagnetic 
properties of the superconducting material are described by a pair ((f), d) which minimizes the free energy 
functional 

g = l J d/ ^ (tf(tf 2) + ^ + 24> 2 d 2 + 2(a' - h e f)dx (i.i) 

The functional Q is now known as the Ginzburg-Landau energy and provides a measure of the difference 
between normal and superconducting states of the material. The variable <j> is the order parameter which 
measures the density of superconducting electrons, and a is the magnetic field potential. Also, h e is the 
external magnetic field, and k is the dimensionless material constant distinguishing different superconductors, 
i.e. < k < -j= for type I superconductors, and k > -j= for type II superconductors (see also [7]). The 
minimizer requirement, that Q be stationary with respect to general first order variations of the functions 4> 
and a, leads to the boundary value problem 

f = fc 2 0(0 2 +5 2 -l), (1.2) 

a" = 4> 2 d, (1.3) 

0'(±|)=O, ~a>(±±)=h e . (1.4) 

It is routine to prove that Q has a smooth minimizer satisfying (1.2) - (1.4) for any positive h e . In 
1964 Marcus [13] investigated the problem (1.2)-(1.4) and gave arguments which imply that a non-trivial 
minimizer of Q should also satisfy 

<f>(x) > for all x£ - ^, ^ , (1.5) 
and therefore this is the only kind of solution we consider. A solution of (1.2)-(1.3) is called symmetric if 

4>'(0) = and 5(0) =0. (1.6) 

It follows from (1.2) and (1.3) that if <j)'{0) = 5(0) = 0, then </> is an even function and a is odd. Thus 
is symmetric with respect to the origin and a is antisymmetric. If (1.6) does not hold then the solution is 
called asymmetric. Marcus makes the conjecture that a minimizer of Q is probably a symmetric solution 
satisfying (1.2)-(1.6). However, he leaves open the possibility that asymmetric solutions may also exist. 

In later work Odeh [15] gave criteria for asymmetric solutions to exist by bifurcation as h e increases, 
and in [2]- [5], Bolley and Helffer give results implying that these criteria are satisfied for each k > 0. The 
existence of at least one symmetric solution has been investigated by Odeh [15], Wang and Yang [18], Yang 
[19], and also Bolley and Helffer ([2]-[5]). Numerical studies, such as the work of Seydel ([16], [17]), and also 
more recent theoretical work of Kwong [11], predict that a range of parameters exists for which asymmetric 
solutions and multiple symmetric solutions coexist. The work of Seydel [16] also predicts that there is a 
range of parameters for which no non-trivial symmetric solutions exist, yet asymmetric solutions do exist. 
Other papers, such as [6], considered the problem on an infinite interval, whereas in our work, the interval is 
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large but finite. None of these studies addressed the physically important criterion of whether the solutions 
of the problem (1.2)-(1.5) are actually global minimizers of the energy functional Q . However, a recent paper 
by Aftalion [1] does discuss this problem, and includes a conjecture that asymmetric solutions can have a 
lower energy than the symmetric solutions. This is confirmed in the current paper. 

In two recent papers ([9], [10]) we began our investigation of the problem (1.2)-(1.5) with the goal of 
proving the existence of solutions predicted by the numerical studies described above. First, we studied the 
existence of multiple symmetric solutions and proved the following result. 

Theorem 1 (see [9]) 

(i) Let k G (0, -j=). If d > is sufficiently large there is a range of values of h e for which at least two 
symmetric solutions exist. 

(ii) Let k > If d > is sufficiently large there is a range of values of h e for which at least three 
symmetric solutions exist. 

Remark: In related studies, Bolley and Helffer ([2]- [5]) have investigated other properties of symmetric 
solutions, including bifurcation analyses and the uniqueness of solutions. Their analysis assumes that k tends 
to zero, whereas our results assume that k > is fixed and d becomes large. 

In [9] we shifted our attention from the symmetric solutions found in [10] to the study of asymmetric 
solutions. As mentioned above, the numerical experiments of Seydel (in particular, see Fig. 6.10 in [16]) 
predict that there is a range of parameters in which there are no non- trivial symmetric solutions, yet asym- 
metric solutions do exist. This leaves open the possibility that in this parameter range the minimizer of Q 
could be an asymmetric solution. Thus, our goal in [9] was to prove that there is a range of parameters in 
which only asymmetric solutions exist, and that the energy Q is minimized by such solutions. The first step 
in proving this result is to find an upper bound on the values of h e for which a symmetric solution exists. 
Thus, for fixed k > and d > 0, we let h* denote the supremum of the set of all positive h e for which a 
non-trivial symmetric solution exists. Since we are assuming that d> 1 wc define h sym = lim^oo^e- 

In [9] we proved that if k > ^.m ' t ^ LCn ^e Vm < V^fc for large d. Unfortunately, a rescaling error led us 



to assert that this inequality was sufficient to prove that there are asymmetric minimizers. It turns out that 
we need a stronger estimate. We shall show below that the inequality 



suffices. This is obtained by a routine but tedious refinement of the proof in [9]. We shall not repeat this 
proof here. However, in an appendix we describe the changes which must be made in [9] to obtain (1.7). We 
add that the proof of (1.7) is considerably easier for large k. The proof for this case is in [9]. However, we 
want to include values in the range of type I superconductors, and this is more difficult. 

Statement of Main Results 

In this paper we will make use of (1.7) in proving that Q has an asymmetric minimizer. We prove 
two main results. First, we fix k > -y== so that both type I and type II superconductors are included. 
Then, in Theorem 2, we consider large ii > 1 and prove that there exists a family of small amplitude 



h s e ym < 1.68k 



(1.7) 
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asymmetric solutions of the problem (1.2)-(1.5). We will also prove that ^ > 1.6831 for large d, for each of 
the asymmetric solutions found in Theorem 2. This and (1.7) confirm the numerical prediction of Seydel in 
[16], that there is a parameter regime in which there are asymmetric solutions, and only trivial symmetric 
solutions. 

The work of Bolley and Hclffcr includes results which imply the existence of asymmetric solutions. The 
proofs, which must be pieced together from several papers, are by bifurcation theory, and do not appear to 
give the estimate of h e which we obtain and which is essential in our discussion of whether these asymmetric 
solutions are energy minimizers. The proof of Theorem 2 is given in Section 2. Next, in Theorem 3 (proved 

in Sections 3 and 4) we show that the energy is negative for the asymmetric solutions found in Theorem 2. 
This will allow us to prove that for large d, a global minimizer of Q must be asymmetric. 

In addition to showing that asymmetric solutions exist and have small amplitude, we will prove that 
each of these has exactly one critical point (a relative maximum) in the open interval (—5, 5), and that the 
relative maximum occurs at a value close to — |. Because of these properties, we find it convenient to rescale 
the problem. We introduce parameters r, h, to and M, and a new independent variable t, by setting 

r=l, h =!±, m + M = kd, ^^Jf^ + rn-M). (1.8) 

Next, we define new dependent variables ip and A by 

4>(x) = (3ip(t) and a{x) = A(t), (1.9) 
where /3 = 0(0). Then (1.1)-(1.5) become 

P 2 f M f,2,0 2 <P 2 



2(m + M) 



j m (^^r 1 + a2) + (v/)2 + w [a ' (lio) 



=iP(/3 2 tP 2 + A 2 -l), (1.11) 
A" = rf3 2 ^ 2 A, (1.12) 
V>'(-m) = V'(M) = 0, A'(-m) = A'(M) = h > 0, (1.13) 

and 

V>>0 in [-m,M]. (1.14) 

We now state our existence result: 

Theorem 2 For sufficiently small (3 > there is a value h > 1.6831 and a solution (ip, A) of (l.ll)-(l.l^) 
defined on an interval [—m,M], such that the following properties hold: 

(i) m = m((3) > 0, M = M(/3) > and lim /3 ^ + (to, M) = (0,oo); 

(ii) ip' > on (-m,0), ip'(0) = and 0(0) = 1; 

(iii) ip' < on (0,M); 

(iv) There is an ft, > 1.6831 such that 
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A'(-m) -^h° as (3 -» 0+ 



(1.15) 



Remark: Because of (1.8) and (1.9), each of the solutions (ip, A) found in Theorem 2 corresponds to a 
solution (<p, a) of (1.2)-(1.5). Since tp has a relative maximum at t = 0, (1.8) implies that ^'(x max ) = 
where x max = f ( ™+m )- It follows from (1.9) and properties (i) and (iii) of Theorem 2 that for small 
(3 > 0, £ max < 0, and < for alia; e (£ max , |). Therefore </>'(0) < 0, and (1.6) cannot hold. We conclude 
that (</>, a) is an asymmetric solution of (1.2)-(1.5). 

We now state our second result. Recall that Q gives the free energy of a solution. 

Theorem 3 Let (tp, A) be an asymmetric solution found in Theorem 2. There exists 7 > such that if 
re (0, 2 + 7) then Q < for sufficiently small (3 > 0. 

Asymmetric Minimizers 

We now return to the original system (1.1)-(1.5) and show that Q has an asymmetric minimizer. Recall 
that k > ^2.01 1S fi XC( I- Also, it follows from (1.7) that if ^ > 1.6831 and d 3> 1 then the only symmetric 
solution of (1.2)-(1.5) is the trivial solution (<fr, d) — (0, h e x), also known as the "normal state." Substitution 
of this pair into (1.1) shows that Q = 0. Next, we conclude from (1.8), (1.9) and the results of Theorems 2 
and 3 that the problem (1.2)-(1.5) has an asymmetric solution for large d, that o'(— f ) > 1.6831 k, and that 
the corresponding energy Q is negative. Therefore, in this parameter range, since Q is zero for the trivial 
symmetric solution and negative for the asymmetric solutions, a minimizer of Q must be asymmetric. 



2. Proof of Theorem 2. 

Our goal in this section is to show that for small [3 > there is a solution (ip, A) of the system 

V>" = V(/3V + ^ 2 -l), (2-la) 

A" = r[3 2 i} 2 A, (2.16) 
on an interval [—to, M], where to > is small and M > is large, and such that 

V>'(-to) = ip'(M) = 0, A'(-m) = A'(M) = h > 0, (2.2a) 



ip 1 > on (-to, 0), V' < on (0, M), (2.2b) 

ip > on [-m,M], (2.2c) 

and 

V>(0) = 1, V'(0) = 0. (2. 2d) 

Our method of proof uses a topological shooting argument. For this we begin by analyzing the important 
properties of solutions of the initial value problem (2.1a), (2.1b), (2. 2d) when [3 = 0. 
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In this case A' is constant, and we set A' = h, where h > is to be determined later. Setting A(0) = Ao, 
also to be determined later, we obtain the second order linear equation 



i> = ((A + htf - (2.3) 
Because of (2. 2d) we consider the solution of (2.3) such that the ip(0) = l,ip'(0) = 0. 

Lemma 2.1 Suppose that —1 < Ao < 0. Then there is a unique ho > (depending continuously on Ao) 
such that ip > 0,ip' < on (0, oo), and ip(t) — > as t — > oo. If < h < ho then ip = before ip' — 0, w/ii/e i/ 
h > ho then ip' = before V = 0. 

Proof: We consider the Ricatti equation obtained by setting 
Then p(0) = and 

where i?(s) = (A + s) 2 — 1. Since A € [— 1,0], p initially decreases, for any h > 0. Further, p' < as 
long as p(s) 2 > H(s). As long as p' < 0, the right side of (2.4) is an increasing and negative function of h 
and p. Suppose, for some first s > 0, that p'(so) = 0. Then H(so) = p(so) 2 so that A(s ) > 1. Therefore, 
p"(so) = 2A< ^°) > 0. For s > s it follows from the equation for p" that p" > so that p increases until 
p = 0; i.e. ip' = 0. Also, if p' becomes positive for some hi, because p crosses the curve p 2 = H(s), then the 
same must happen for any h > hi. To see that there are values of h such that p' becomes positive we refer 
to (2.3). From that equation and our assumption that ip(0) — 1, ^'(0) = we easily see that for large h, ip' 
becomes positive before ip = 0. This implies that p, and hence p' , must become positive for large h. It then 
follows from continuity that the set of values of h > such that this happens is open. 

For small ft, > 0, on the other hand, we see that p decreases to below the curve p — —s. For example, if 
— 1 < Ao < and < h < 1 — Aq , then p'(0) < — 1, so immediately p decreases below — s. It Ao = — 1, then 
p(0) = p'(0) = 0. However, p"(s) < ^j^- and this integrates to show that for small h, p(2) < —2 and p' < 
over (0,2]. It follows from (2.4), and our assumption that A £ [—1,0], that if p' continues to decrease until 
p — > — oo at some finite si. That is, ip(si) = 0. It is clear that the set of h's such that p(s) < — s for some 
s is an open subset of the interval < h < oo. Similarly, as we pointed out above, the set of h's such that 
p(s) 2 < H(s) for some s is also open. 

Moreover, if p ever falls below — s, then we see from the equation for p" that thereafter, p' < —1, and 
p decreases monotonically to — oo. On the other had, if p' is ever positive, then we saw that p becomes 
positive, and because of (2.4) it must remain positive. 

Hence there is at least one ho such that on < s < oo, p'<0 and > p > — s. These bounds imply that 
the solution exists on the entire interval [0, oo). Further properties of this solution are given in the following 
result. 



Lemma 2.2 For such an ho, p' < 0, — s < p(s) on (0, oo), p(s) < —^/H(s) if H(s) > 0, and p + ^/H(s) — > 
as s — > oo. 
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Proof: The inequalities have already been proved. To see the limiting behavior of p(s), we first note that 
for h — ho,p'(s) is bounded. If p' is unbounded, then it must get arbitrarily large and negative, but then 
the equation for p shows that p' remains large and negative and p(s) < — s for some s. 

Writing p' = W2LzgM^2I±g) ana ; noting that the first of these factors is unbounded, we see that y/H(s) + 
p(s) tends to zero, which implies Lemma 2.2. 

To prove uniqueness we suppose that there is a second positive value of ft, hi < fto for which this behavior 
occurs. Let pi, po denote the corresponding solutions. Then pi(0) = po(0) = 0, and it is easily shown that 

^(Pi-Po)<0 

for all s > 0. Thus pi — p could not approach zero as s — > oo. However, Lemma 2.2 implies that p\ — p = 
(pi + VH) — (po + VH) — > as s — > oo, a contradiction. Therefore ho is unique. 

To complete the proof of Lemma 2.1 we must show that ho depends continuously on A n . This follows 
from the uniqueness of h by a standard argument. 

Now let A = —1, and define 

ft = lim ho(0). 

We will find the desired solutions near the point (A , h) = (—1, ft ). 
Lemma 2.3 ft > 1.6831. 

Proof:Let m m 

ptt) = w 

so that 

p' = H(ht) - p 2 . (2.5) 

Here, because A a = —1, H{s) = —2s + s 2 , where s = ht as before. We use the following result about a 
solution p of (2.5) such that p(0) — 0. If p " < for some ti, with p,p',p" < on (0, t\) then p decreases 
monotonically to — oo. To see this, compute the equation satisfied by p . That is, 

//// /// // 

p + 2/5/5 = -6/5 p . 

If there were a first y > ti where p (y) = then p (y) > 0. However, the last equation gives p (y) < 
since the definition of y implies that p'p is positive and increasing on (t\ , y) . From these derivative properties 
it follows that p decreases below — s, and hence p(s)— > — oo at a finite value of s. Next, we define a sequence 
{Pn} of functions, all defined on [0, 2] as follows: 

p {t)= f {h 2 r 2 -2hr)dr , p N+1 (t) = p (t) - f p N {rfdr for N > 1. 
Jo Jo 

From this definition it is evident that {pn} forms a decreasing sequence of functions defined on [0,2], for all 
N > 1. Furthermore, it is easily shown that our solution p(t) satisfies p < p i on [0,2], for all % > 1. Setting 
ft = iqqq ; we used the computer algebra program Maple to compute p§, a polynomial of degree 225. All 
calculations are with integers and rational numbers so that there are no roundoff errors. From pg and its 
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first two derivatives we can compute p'", without having to compute p-j. We find that p"'(6/5) < 0, and that 
on [0,6/5], p,p',p are all negative. This completes the proof of the lemma. 

Our solution ("4>,A) to (2.1a)-(2.2d) is obtained by perturbing A(0) from —1, keeping (3 — 0, and then 
letting (3 be positive. Our argument is by "shooting" , rather than by use of bifurcation theory. 
Thus, we will assume that 

A(0) = -l + e 

for small e > 0, and we let ho = ho(e) denote the value of h found in Lemma 2.1. 

In constructing solutions to (2.1)-(2.2), in order to get something meaningful at (3 — 0, we replace the 
boundary conditions on A in (2.2a) with 

/M 
iP 2 (t)A(t)dt = 0. 
-m 

Lemma 2.4 Suppose that e = 0, (i.e. Aq = —1) and let tpo be the solution found in Lemma 2.1 where 
h = h° = h (0). Then 

/ iP (t) 2 (-l + h°t)dt = Q. 
Jo 



Proof: It is easily seen that ipo(t) — > exponentially fast as t — > oo. Thus, the integral in the lemma 
converges and lim t ^ 00 (— 1 + h°t) 2 ipo(t) 2 = 0. We set ip = tpo in (2.3), multiply by ip' , and integrate by parts 
to obtain the result. 

From Lemma 2.1 and the definition of ho(e) we see that for each h > ho(e), there is a first t = th > 
such that tp'(th) — 0,tp' < and ip > on (0,th), with ip(th) > 0. By the implicit function theorem, th is 
continuous in h on (ho(e), oo) since ip (th) > 0. (If ip (th) = then th = \ and ip (th) > 0, a contradiction.) 
Further, th — > oo as h — > h (e) from above. 

Now consider small e > 0, set h = h (e), and compute ip" (0) and ^"'(0). We find that ip"(0) = — 2e + e 2 
and V"'(0) = 2(— 1 + e)feo( e )- From this and the fact that /io(e) — > ^° > as e — > + , we see that for small 
e, tp'(—m) = for some m = m (e) > 0. Furthermore, 

-m= T ^ + ^ 2 )- (2.6) 
Here, 0(e 2 ) < Le 2 for some L independent of e. We now compute 

/oo 
-m 

where A(s) = — 1 + e + ho(e)s. Again, multiplying (2.3) by ip f and integrating by parts gives 

/oo 
iP(s) 2 A(s)ds. 
-m 

Using (2.6) and the fact that ip(-m) = 1 + O(e), wc find that 1(e) = -f- + 0(e 2 ) as e -> 0+. Thus, we 
have shown that with h = ho(e) for small e, /(e) < 0. Fix e > small enough that this inequality holds 
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for h — ho(e). Then for h — h (e) positive but sufficiently small, — m = —m(e, h) will still be defined as the 
largest negative zero of ip', A is positive on (th, oo), and therefore 

[ h ip 2 Ads < 0. 

J —m 

Our goal now is to find an h and e such that — m and t h are defined as the negative and positive zeros 
of ip' closest to t = 0, with m small and th large, and such that 

I=( %j} 2 Ads > 0. 

J —m 

We will see that it is not necessary to have the same m and th as before. Starting with e = and h = ho(e), 
we now raise h, instead of e. Thus, for small h — h a > our solution satisfies A(0) = — 1, A'(0) = h > h = 
ho(0), ip'(th) = for some large th- In this case, we multiply (2.3) by ip' and integrate from to th, where 
th is the first positive zero of tp' . Using integration by parts once again, we find that I > 0. We then keep h 
fixed and raise e slightly, whereupon both m = m(e, h) and th are defined, still with / > 0. We summarize 
these results in the following lemma. 

Lemma 2.5 For each sufficiently small e > there is an h\{e) > h n (e) such that if ho(e) < h < h\{e), then 
the solution of (2.3) with A n = — 1 + e, ip(0) = 1, and ip'(0) = is decreasing on an interval [0, M], increasing 
on an interval [— m, 0], and ip'(-m) = tp'(M) = 0, and 

r-M 

I(e, h)= tp 2 Ads < 0. 

J —m 

Furthermore, ip > on [— m, M] and hi(e) — > ft-o(O) = h as e — > + . On the other hand, for each h > ho(0) 
sufficiently close to ho(0), there is an interval (0, ei(/i)) of e's such that ip has the same behavior, but 
I(e, h) > 0. As e and h — h (0) tend to zero, m — > and M — > oo. 

Corollary 2.6 For e > sufficiently small, there is an hi = /12(e) sucA t/iat ifip is the solution of (2.3) with 
A = — 1 + e and (ip{0), ip'{0)) = (1,0), i/ien t/iere are values M > m > suc/i t/iat 

V/(-m) = tp'(M) = 0, V' > on (-m, 0), -0' < on (0, M), (2.7a) 
V>"(-m) 7^ 0,V"(M) ^ 0. (2.76) 
Further, I(e, /11(e)) < and J(e, /12(e)) > 0. Finally, /12(e) — ► /io(0) = /i° as e — > 0. 

With e > sufficiently small so that Lemma 2.5 and Corollary 2.6 hold, we now raise (3 and consider 
solutions of (2.1) satisfying (2. 2d). Since ip" 7^ at the zeros of ip' when (3 = 0, it follows from the Implicit 
Function Theorem that the zeros — m < and M > of ip' persist as continuous functions of [3 and the 
values of A(0) and ^4/(0). For each sufficiently small e > there exist functions /11(e) and /12(e) independent 
of (3, such that for sufficiently small (3 (depending on e) the solution (ip,A) of (2.1) and (2. 2d) with 

A(0) = -l + e, A'{0) = h 1 (e), (2.8) 
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satisfies (2.7) on an interval [—m,M], and 



/ 



ip z Adt < 0. 



Also, the solution of (2.1), (2. 2d) with A(0) = -1 + e, A'(0) = h 2 {e) satisfies 

i-M 



ip 2 Adt > 0. 

J —m 

Furthermore, — > as e — » 0, and 

(m,M) - (0,oo), (/ii(e),/ l2 (e))^(/ i °,/i ). 

In addition, it follows from (2.1b) that A is uniformly bounded on [— m, 0], and therefore m) = h — > ft 
as e — > and /3 — > + . It then follows from continuity that for given fixed small e, and sufficiently small 
> 0, there is an /i e (/ii(e), /12(e)) such that the solution (ip e ^,A €i p) of (2.1), (2. 2d), and (2.8) with 
A'(— m) = h also satisfies J M m ip 2 Adt = 0, so A'(M) = h. The conclusion of Theorem 2 now follows from 
the transformation (1.8)-(1.9). 

In the next section we turn to the proof of Theorem 3. For this we will use the following result. Let 



F(ip,A)(t)dt 

-m 



where 



a M N 2 

^( S fA{ s )ds) - 



Lemma 2.7. With (ip, A) = (ip €t p,A €t p) chosen as in the proof of Theorem 2, we have 

\im{Yim£^,A)}=£^ ,A°) 

where A°(s) = -1 + h Q s. 

Proof: Our proof of Theorem 2 shows that 

lim{Um(^( S ),^( S ))} = (Ms),A°(s)) 
uniformly on compact intervals. Let S > be given. We can choose K\ > such that 

ti>(K{f<5 / \F(1> ,A°)(t)\dt<5, (2.9) 

and 

A°(s) 2 -2> (2.10) 
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for s > K\. Then for sufficiently small e > we choose j3\ — /3i(e) > such that for < f3 < f3\, 

[ Kl (F(^A)(t) - Ftyo,A°)(t))dt 
Jo 



< S. (2.11) 



Further, we can insure that for < (3 < (3\, (2.10) also holds on [K\, M] with A substituted for A . 
We consider two cases: 

ffill > -1 (i) 

and 

Let p = so that, from (2.1a), 

As) > Ms? -I-P 2 - 

In case (i) we have p' > A ^ on [Ki, M] as long as — 1 < p < 0, so that p = (and hence ip' = 0) before 



/ — f/.s • = 1. 

In other words, 



> Kl 2 



M Af \2 



A{S) ds < 1. (2.12) 



Ki 2 

Thus, (2.10) implies that if M/get > K x , then ^-2)/ gel on [t,M], and hence 



J iP(s) 2 A(s)ds < V(^i) 2 J -f-dsK^K^ 



it Jt 
Since (2.10) and (2.12) also imply that M - K x < 1, it follows from (2.12) and (2.9) that in case (i), 

i-M 

/ F(ip,A)(t)dt < 5 (2.13) 

JKi 

for sufficiently small e and (3. In case (ii), since ip'(M) = 0, there must be a T e [K\,M) such that p < — 1 
on [ifi,T] and p > -1 on [T, M] . First consider F(ip, A)(t)dt. On [ifi,T] we have V' < -V>, so that 

^(t) < V(^i)e Kl ~*. (2.14) 

We now estimate the term 



/LVi 
iP{s) 2 A(s)ds 

in F(ip, A) for K x < t < T. For this we use: 
Lemma 2.9 lim^^ ip (i)A°(t) = 0. 



Proof: In our estimates of ipo we saw that + flo(t) -» 0. Since A grows linearly, the result follows. 
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We now continue with the proof of Lemma 2.7. We choose K\ so that in addition to the earlier constraints, 
ipo(Ki)A°(Ki) < 1. For small e and (3, this inequality will also be satisfied by (tp, A). If t < T, then 

/M pT pM 

ip(s) 2 A(s)ds = J ^{s) 2 A{s)ds + J iP(s) 2 A(s)ds 

pT pM 

< / iP{sfA{s)ds + tp(Tf / A(s)ds. 



The argument used earlier to get (2.12) shows that J^f A(s)ds < 1, where we use (2.10) with A substituted 
for j4°. From an integration of ^ < -1 we obtain 

/ iP(s) 2 A(s)ds < tp(t)( max ip(s)A(s)) [ e^ s ds. 

Jt t<s<T J t 

In [K U T], ip' < -ip, so £{i>(s)A(s))< -ip(s)A(s) +ip(s)A'(s). But A'(s) < A'{M) = h < 2h° for small e 
and P, so ip(s)A(s) is decreasing in [K\,T]. Hence, ip(s) 2 A(s)ds < ip(t) < tl)(K\)e Kl ^ t . Using this we 
find that on [K U T], 

F{^,A) < r^(^i) 2 e 2(Kl - 4) . 

Thus, we can further restrict K\ to insure that for 
small positive e and /?, 

/ F(il>,A)(t)dt<5. 
Finally, we consider f£? F(ip, A)(t)dt. As in case (i), we have 



pM 
JT 



2 

Hence wc find that 



A(s) 2 , 

; ds<l. 



/M 
i(;(s) 2 A(s)ds < ^{T f < ^(K^ 2 



and 

M-T < 1. 

We then have, (without further change in K{) that 

pM 



pm 

J F(i/),A)(t)dt<6. 



IT 

Since S was arbitrary, this completes the proof of Lemma 2.7. 
3. Proof of Theorem 3 

In order to prove Theorem 3 we need to show that the asymmetric solutions found in Theorem 2 have 
the additional property that their corresponding energy Q is negative if e > and j3 > are small enough. 
Thus, wc define 
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Q = ?^ S . (3-D 

Then, by (3.1) and (1.10), Q is given by 

Q = ^J M m (f(rf - l + A2 ) + ^ + h & L ( 3 - 2 ) 

Therefore, if we show that Q < for small e > and (i > then Q also is negative and Theorem 3 is proved. 
We need to simplify Q. For this we use (1.11) and conclude that 

^(^f 1 + A 2 ) + m 2 = + - ^f- (3-3) 
Then substitution of (3.3) ino (3.2), together with the observation that (tpip'Y = ip^) + {ip') 2 , reduces Q to 

1 r M B 2 ib 4 1 

Q = J JW)' -^r + ^ A '- h ^ dt ^ 

Since ip'(-m) = ip'(M) = 0, (3.4) further reduces to 
Next, it follows from an integration of (1.12) that 

/M 
ifAds. (3.6) 

Finally, we substitute (3.6) into (3.5) and arrive at 

(r(y tfAds) -\)dt. (3.7) 

In view of Lemma 2.7 we conclude from (3.7) that 

f°° / / f°° \ 2 tl) A \ 

*3>]& Q = J Hi <* A da ) " 2)* (3 ' 8) 

where A = — 1 + h°s. 

In the next section we prove 

Lemma 3.1. There is a value 7 > such that 

r°° / r°° \ 2 ,/,4 

jf (r(jf Vo 2 ^°^J -f)dt<0 (3.9) 

for all r £ (0,2 + 7). 

From Lemma 3.1 we observe that if r £ (0, 2 + 7) then Q < for small e > and /3 > 0. Therefore Q is 
also negative for small (3 > and Theorem 3 is proved. 
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4. Proof of Lemma 3.1 

For the proof of Lemma 3.1 we recall that Vo satisfies 

%=M-2h°t+(h°) 2 t 2 ), (4.1) 

Vo(0) = 1,^(0) = (4.2) 
where h° > is the unique positive value for which tp satisfies ip' < for all t > 0, and 

lim(^ (i),VoW) = (0,0). (4.3) 

The existence and uniqueness of h° was proved in Lemma 2.1. Thus, our main objective is to prove 
Lemma 4.1 There is a value 7 > such that 

(r(J% 2 (-l + h° S )d s y -^(t))dt<0 (4.4) 

for all r G (0,2 + 7). 

The proof of Lemma 4.1 relics on an auxiliary result which we now establish. 
Define the Ricatti variable q — ^jfy- Then q satisfies q(0) = and 

q' + q 2 + 2h°t - (h°) 2 t 2 = 0. (4.5) 

g(0) = q'(0) = (4.6) 
Lemma 4.2 It follows from (4-5) and the definition of q that 

q<0 and q' < for all t > 0, (4.7) 



/ 

Jo 



q < -^(h°)H 2 -2h°t for all t > ^, (4.8) 

q' > -h° for all t > 0. (4.9) 

Proof: Setting h = h° and s = th°, we observe that (4.7) and (4.8) follow immediately from Lemma 2.2. 
It remains to prove (4.9). It follows from (4.6) that ^'(0) = 0. Thus q' > —h a on an interval [0,77) for small 
n > 0. Suppose that (4.9) is false. Then there is a first t > for which q'(i) = —h°, and q (i) < 0. Two 
differentiations of (4.5) lead to 

q" +2qq" =2((h°) 2 -(q') 2 )- (4.10) 

One solution of (4.10) is q' = —h° for all t. Thus, if q (i) — then uniqueness of solutions implies that 
q' = —h° for all t > 0, and in particular, q'(0) = —h°, contradicting the fact that q'(0) = 0. Therefore it 
must be the case that 

q'{i)=h° and q" (t) < 0. (4.11) 



14 



It then follows from (4.10) and (4.11) that q" (t) < for all t > t so that 

lim q'(t) < -h°. (4.12) 

We conclude from (4.12) that 

lim q(t) + h°t < 0. (4.13) 

t — >oo 

However, it follows from Lemma 2.2 that lim^,^ q(t) + h a t = 0, contradicting (4.13). Thus, it must be the 
case that q' > —h° for all t > and the proof is complete. 

We now return to the proof of Lemma 4.1. First, we conclude from (4.7) and (4.8), and the properties 
that < Vo < 1 and ip' < for all t > 0, that the integral ipQ(t)dt is well defined and positive. The 
same reasoning shows that 

J = J (J ^(s)(-l + h°s)dsj dt (4.14) 
is well defined and positive. It remains to prove (4.4). For this we begin by estimating the integral 

/oo 
ipl(-l + h°s)ds. (4.15) 

An integration by parts reduces (4.15) to 

?/; 2 1 f°° 

H = + h°t) 2 - Yo J t V^o(-l + h^fds (4.16) 

It follows from (4.1) that V'oV'o + V'oV'o = V'oV'o( — 1 + h°t) 2 and therefore (4.16) becomes 

"=-i(-i+*°<) 2 + ffi + 5>;f. <«-) 

Recall that V' = qip- Then (4.6) and (4.17) imply that 

It follows from (4.18) and (4.14) that 

J = £ H\t)dt< ^L{q'Ydt. (4.19) 
Because Lemma 4.2 gives — h° < q' < for all t > 0, (4.19) further reduces to 

J< lj o ^(Wt. (4.20) 
Finally, it follows from (4.4), (4.14) and (4.20) that there exists 7 > such that 

poo poo 2 /.4 * 

Jo v(i ^o(-^ + h°s)ds) -f)dt<0 
for all r £ [0, 2 + 7). This completes the proof of Lemma 4.1. 
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Appendix: In [9], when studying symmetric solutions, we used a different scaling from that used in this 
paper. That is, we rescaled (1.2)-(1.5) by setting 

T ^ k 2 d 2 , h e d 1 2x 

and defined new dependent variables (f> and a by 

4>{y) = 4>(x),a(y) = a(x). 
This transforms the problem (1.2)-(1.5) into 

cj>" = K<t>{tf +a 2 -I), 
a" = rK(f> 2 a, 

<p > on [-1,1]. 

Since we are considering symmetric solutions, we consider the initial values 

0(0) = 0, <t>'(0) = 0, a(0) = 0, a'(0) = a. 

In [11], Kwong proved that for each [3 G (0, 1) there exists a unique a = a{j3) > 0, continuously dependent 
and decreasing in f3, such that <f>' (1, f3,a((3)) = 0. He then set h(/3) = a'(l, f3, a(/3)) and proved that h is 
continuous, h(0) > and h(l) = 0. Thus, to obtain the upper bound h sym < \/3 given in [9] we found that 
it was sufficient to fix r e (0,2.01] and estimate 

h(py 



h sym = limK^oo ( sup —7=). 

V 0</3<1 vK ' 



In [9], our estimate for h sym was obtained by carefully analyzing the behavior of (<j>, a) for each f3 € (0, 1). 
We considered three intervals of (3, namely (0, .1], (.1, 1 — -^=\ and (1 — 1). We made repeated use of the 
"energy" function 



A/2 a ,2 



2 Ji2 



1 V <f> - a 

K rK 2 v 

which is constant for solutions of the system. We defined the function 

and found the point yo where Q(yo) — 1 to be important. The difficult part of the estimate required us to 

consider a small interval of values of 4>(yo), say I\ < 4>(yo) < ^2- A priori we know only that < (j)(yo) < 1, 
and to get better estimates we had to subdivide (0, 1) into nine small subintervals [Ji, I2]. 
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In this paper, to get the improved estimates required to show that h sym < 1.68 simply requires that 
we use more subintervals. Thirty-four subintervals suffice. They are: [0,.407], [.407, .56], [.56, .6], then thirty 
intervals from .6 to .9 in steps of .01, and finally [.9,1.0]. 

With this change, the proof in [9] gives the required upper bound for h sym , and because of our rescaling 
this immediately leads to hl ym < 1.68fc for k > ^ m ■ As we mentioned in [9], our estimate for h sym is much 
easier for small r. 
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